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NEUTRON DIFFUSION AND RESONANCE EFFECTS ON FISSIONABLE 
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By L.P.Abagyan , N.OMazazyantz, UJIondarenko , 0. A .Guseynov, A. A. Lukyanov, V ht.Mahanov V.lMe- 
lentyev , U.N .Nikolaev, V .V .Orlov, N .N.Rabotnov, A.P.Souvorov , L.N.Usatchyov , V. V. Filippov. * 

The importance of taking into account the resonance structure of nuclear cross sections 
(especially for U^8) was clearly understood in early investigations on thermal fission 
theory [1-6]. 1. 1. Bondarenko pointed out in 1957 that resonance effectsare of importance 
in the fast reactor physics too. Later it was confirmed by macroscopic experiments and cor- 
responding calculations. 

In this report the general results of the resonance effects investigations and methods of 
its applications to the purposes of reactor calculations are presented. Some of these calcula- 
tions have been already reported [7-15] and used in the construction of a multigroup constant 
set. The theoretical and experimental problems of resonance capture in heterogeneous thermal 
reactors are not discussed. 

RESONANCE STRUCTURE OF NUCLEAR CROSS SECTIONS AND 
REACTOR CALCULATIONS 

When considering the energy spectrum of neutrons diffusing in the medii two points could 
be marked. 

First, the neutron flux varies sharply reflecting resonance structure of the cross sections; 
these variations form the ’’Microstructure” of the spectrum. Second, the neutron absorption, ' 
diffusion and existence of neutron sources result in rather smooth changing of average neutron 
flux and its deflections from Fermi spectrum. A number of methods could be used for the 
determination of neutron space-energy distributions in the case of medii with smoothly varying 
cross sections. As for the calculation of neutron spectrum ’’microstructure”, the exact solu- 
tion of the problem depends on the detailed information about the cross-section variations in 
the wide range of neutron energies and thus seems to be impossible in present. 

Fortunately, the purpose of calculations usually is to determine the values averaged 
over the great number of resonance. (Sometimes the exoct information about the neutron flux 
is nevertheless required. But a number of such rosonar.es is not too large and an individual 
consideration does not make a trouble). For this aim the knowledge of average neutron flux 
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0 (?, u, n) is sufficient. The evnluntion of the number of events corresponding to uny 
process can be carried out by using the cross-section values averaged with neutron flux 
as a weight. 

2 X = [ /2 x (u)<£(r, U, sfydu]// <^(r , u,fl)ilu (1) 

Au 

(Index M x” denotes the type of interaction - capture, fission etc.). 

These average cross sections are in general dependent on the direction and coordinates^ 
and their evaluation requires the information about the dependence of neutron flux <^(r,u,n), 
which is in general available. If the diffusion and (or) energy interval used for averaging 
are so small that 

lD A 2<£/<jJl «2 a + £V Au (2) 

(but Au is at the same time great enough for average cross sections to be smooth) the flux 
(£(r, u, O) in the case of one-dimensional geometry has a form 

</>(? u , 11) = (l/4rr) 2(21 + X) (u,? ) </>j(u)pj(/i) (3) 

Here ^ is smooth function U and 0j(u) — (/>j(u, r) /<jT j(u,r) describes the microstructure of 
spherical harmonic spectra and depends on the medium's properties only. Notice that 
although relation (2) is the condition of the applicability of diffusion approximation the 
results are valid also for higher approximations. The condition (2) usually holds if regions 
far from the boundaries as well as from the localized neutron sources are under considera- 
te* >»ecan use available average constants for the description of neutron propagation in 
.egioir i For more accurate determination of the neutron flux distribution near the 
boundaries special consideration is needed. It is true also for the problem of neutron .flux 
determination in theca 3 e of medii with deep interference minima in total cross section when 
condition (2) also fails to hold in principal. We restrict ourselves to the region of problems 
(which is wide enough) where condition (2) holds. 

KINETIC EQUATION IN 1 - REPRESENTATION 
The system of equation for moments (x, u) has a form (for simplicity the plane geometry 
only is considered) 

<t, 0 l + d = Q(x lU ) + If 9 i 0 ^ i 0 (u 1 -* u )du 1 + /^ o X in0 ( u 1 -u)du 1 

0 e 2 + 1(1 + Dd^/dx+ld^/dx I = 2/ (uUu)dul + /^Si^u^ujdu 1 (4) 

Here 1 = 1, 2, , 1 - is the total cross section, 2^, is the i — th component of 

medium and inelastic cross section. As we supposed the quantity has a little value 

and the integrals in (4) as well as quantity Q are the smooth energy functions. That is why 
we may adopt the next expression for neutron flux': <£ 0 = T/2(u)^l/X(u) 
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. V 1 

9o = prjry- where H* - collision density. 


In the same way 


*l<u)-£ a. /S)"'*" 1 (u) 
n=0 1 


where coefficients are determined by recursion formulae': 


1-1 




Herc 2 a ^-S/2 a ,(uU)du 

By averaging liqs (4) we have 


(5) 


( 6 ) 


+ + Sin^u^ujdu 1 


^1 “ 1 + 2T7r ,(1 + 1)(5 ^l + l /ax) + 1( ^-l /<?x)l = p<6 1 2 sl (u 1 -u)du 1 + /^ I 2 ino ( u l-.u)du 1 

(7) 

where averaged cross sections are determined by the relations 

f= <^2'>/<^ l 1 >, 2^ e (u I -u)= <^(u 1 -.u)0 1 (u 1 )>/<,i 1 (u 1 )> 


( 8 ) 


2 x = <2 x /2 > /<l/S > , 5! inl (u 1 -u) = <2 in (u 1 -u) ^ (u 1 ) ■> /<0j ( u 1 )> 


and (u) is determined by Eq. (5). 

The discussed averaging technique is based upon the assumptions (in addition to 

1. The collision density variations in averaging interval Au hecause of Plachek spectrum 
oscillations following every resonance are ignored. These oscillations change microstmcture 
of neutron spectrum (in particular 0 o (u) differs from 1/2 (u) which can be taken into account 
without difficulties by discussed method. But this method is hardly important for the range of 
neutron energies where statistical description of resonance structure is valid. Really the 
contribution of every single resonance into the collision density is comparatively small. 

By taking into consideration the level spacing distribution and rather sharp variation with 
energy this spacing divided by average energy loss one gets more smooth collision density 

This effect may he significant only for a few low lying resonances which require individual 
treatment. 
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2. Resonances are supposed to be narrow relative to the scattered neutron energy loss. 

As a rule this condition is invalid only for a few first resonances which are to be considered 
individually. 

3. For the determination of spectrum (j>\ (u) it had Au been supposed in Eqs (4) 

(uLuJdu 1 ~<0,X s /I 

which could be not true in the case of* large energy losses (inelastic scattering and scattering 
by hydrogen). By using iteration method one could calculate the spectrum $ j (u) more accurately. 

GRILLING - HERTZBL APPROXIMATION 

There are systems where an anisotropy of neutron flux is appreciable only for high 
neutron energies. Tn this energy range the inelastic scattering contributes mainly into modera- 
tion. In such a situation the elastic moderation could be determined by using the Grilling- 
Hertzel approximation leading to the next equations 

a /*» I (1-a y)/du) 2 + yS oif ] ^ - d^/dx + Q / 0 o S in .f ino ]l + -* 0 S aif = Q - d<f>/dx + (<t> 0 l in f iao 

V 1 ] + 1/3 i«?0 o /*c+2£> 0 2 /axi = 

^1 [2 1 ■ VP 7T ’ 1 + 1 +(<*) + (1 + lKa^.j/ax!-/ ^.j^du 1 

1 = 2,3 

The elastic moderation in the case of light nuclei can be treated in the same way as the 
inelastic scattering. 

Denote the scattering function /(/Zq , u, u^-»u) moments as 

x en ( u ) - /(u^ •» u) n du^/Pj(g t u , u^- u)dn (10) 

Then by using Eqs (8) the coefficients of Eq (9) can be written in the next form 

£ - < 2 xli 2 • /£'>/ <2 /£■>, y =< 2 x n9 ^ s i ,:> /2 > / < 2 s / S > 

1 UI si i u “ 

J r<piiVi > / < Vf' f/1 =<Vi4 x L>/<V 

Moments x en can be found [II) from the next formulae*: for hydrogen x a =l, x^q^IS, x^=419, 

x 02 = ^ ^ or nuc ^ e * A » 1 
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x 01 = 2/A[l-x in (u)], x„ 9 (u) = 16/3A“)[1 - 213x in (u) + l/2x on (u)] 


10' 


02 


10 v 


20 ' 


* n (u) = <-2/3A)[l— 3x 10 (ul + 2 x 20 (u)] 


( 12 ) 


Where x 10 anc ^ *90 are inoments of angular distribution which is to be found in neutron 
scattering experiment. 

As the averaged in this way cross sections and neutron flux vary smoothly with energy any 
finite-difference method such as multigroup one may be used. These techniques are developed 
for the case of smooth cross sections [17]. We will not consider such methods. 


BOLTZMAN KINETIC EQUATION 

In addition to the solution of transport problem in the spherical harmonic representation 
methods of direct numerical solution of Boltzman Kinetic Equation have been developed 
(S n - method [27], Vladimirov's method [28]. Consider such an equation for neutron flux 
averaged over the resonances. For this purpose Eqs (7) may be used. Multiplying Eq. 1 by 
P}(/i)(2l+ l)/4rr and summingall the quantities we have 

1? (J2I + 1)/ 4rr] Pj (n) + n(d/dx) <j) = Q+S.J\^ (x,u^,/x^)^ s .(u^-+u, /i Q )du^dfi (x,u\g) 2* n (u^ . 


-U, dfl 1 


Where 


f .<»! 


s y ^ =0 (2UD/ 4 ^ s i i(^u)p i ( (Jo 


aR d S. n in the same way 

The sum in the left side of the latter equation may be transformed in the next manner. 
Let I] -> 2 if 1 -* oo. Then this sum is 


— 0 OO 1=0 £ oo 

(14) 

S 1 ^) = 2 [(2l + A)/4ir] ( 2, — 2)Pi((i 0 ) 

1=0 oo 


Where 


Thus the effect of the total cross-section resonance structure is shown to be equivalent 
to some additional neutron elastic scattering described by cross section [-X^ (u f /i 0 )]. The 
total cross section is to be renormalized simultaneously. 

NOtiCe tllat S„ = lim Uoo^l - lim U„/<<M>/<^> 

For example in the case of isotropic elastic scattering by l/2* + * and 

2 = lira, <1/2 1 >/<1/2 1 + 1 > = 2 . 
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where 1 . - the minimum value of cross section in the interval of averaging.- 
min 

THE AVERAGING OF MICROSCOPIC CROSS SECTIONS FOR 
SOME ELEMENTS 

When tne great number of reactor system calculations is required it is desirable to have 
the set of averaged cross-sections for elements used usually in reactor construction. Then 
macroscopic cross sections for any real reactor system can be easily found. Thus the question 
arises: in what manner the cross sections are to be averaged for giving the correct macroscopic 
characteristics of any combination of elements with resonance cross-section structure. This 
problem can be solved exactly if only one of the mixture components shows resonance cross 
section structure within any interval of the averaging. In this case one could find the correct 
average macroscopic cross section Eqs (8) if the values <l/£ n >, <^ s i( u ^ u)/2 n >/<2 n -> 
are found in the next way: <l/2%"^^ n - 2Ni<a*>n 

<S s |( ul -* «)/2 n >/<l/2^> =?N.<4 ) (u 1 -u)>n (15) 

where 

<oS n = <1/(7- +ff) n >"^ n — a % <or si (u* ^ u )> n = <or sl(u^ “*u)/(^ + cr) n >/l/<cA + o)^> (16) 

N- is nuclear density of M i M -th mixture component, a — is total cross section of the other 
mixture components calculated per one nucleus of the considered type. 

When there are more then one element in mixture with resonance cross section structure the 
values averaged macroscopic cross sections depend on the location of resonances in every 
element. In this situation only approximate representation of averaged mixture cross sections 
can be found in terms of the averaged cross sections of particular elements. The uncertainty 
of such a representation would be the least one if when averaging macroscopic cross sections 
(Eqs 15) the average total cross section for all the other elements calculated per nucleus of 
considered element is taken as a. The practice shows u\ ibis case the uncertainty is usually 
less then that of the experimental data. 

The tabulation of average cross section of individual elements for any constant "deludon* 1 
cross section helps a the calculation of macroscopic characteristics for any mixture of such 
elements [16]. 

THE SUBGROUP METHOD 

The discussed above method of takinginto account the resonance cross section structure 
based on the assumption that neutron structure remains constant throughout the medium to be 
considered. Inspite of the approximation is shown to be proved in a number of cases there are 
such a situation when a neutron spectrum structure as well as averaged cross sections Eqs (9) 
are dependent on coordinates. Such a situation always takes place if the contribution of neutrons 
which undergo no collisions (sometimes only some collisions) is not negligible. (For example, 
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the regions near the neutron sources or the medium boundary; the energies corresponding to 
the deep minima in the total cross section curve). Then the information about the cross 
section structure which is contained in average cross sections dependence (Fqs ( 12 ) on the 
’’delution" cross section becomes unsufficient. For the introduction of more complete informa- 
tion required for aforementioned problems such a called "subgroup method" [14] could be 
used. Suppose the smooth energy dependence of neutron flux is described by multigroup method. 
Consider the neutron group corresponding energy interval Ali. 

Let us divide these neutrons into some "subgroups" unifying within any subgroup the 
neutrons with the total cross section value enclosed in small intervals Aojj. Then assume all 
the neutrons of the considered subgroup having the same - average - total cross section as 
well as differential scattering cross section and reaction cross section. If the intervals (Ao^) 
are small enough, the assumed total cross section will be rather close to real one. The same 
situation to some degree will take place for the partial cross sections because of the correla- 
tion with total cross section. The propagation of considered group in continuous medium with 
small gradient can be described by using conventional methods with the help of the constants 
(Fqs ( 12 )) calculated by averaging within the subgroups. 

If the neutron spectrum structure is not constant (and so the subgroup contribution depends 
on the coordinates), one have to calculate the propagation of every subgroup apart. This 
requires the additional knowledge of group-to-group transition cross sections because of the 
elastic scattering, inelastic scattering and fission. As for fission and inelastic scattering the 
problem does not make any difficulty because the resulting energy changing is much more 
than level spacing. But in the case of elastic scattering transition one have to take into 
consideration the relation between the energy loss caused by scattering and level spacing. 

In the isolated resonances region the transition cross sections could be obtained from the 
detailed cross section structure. For higher energies the statistical description based on 
average resonance parameters must be used. When the energy loss caused by scattering 
exceeds the level spacing the elastic transition probability becomes proportional to the contri- 
butions of those subgroups into which the transition takes place. 

Neutron transitions between all subgroups of given group complicate calculations consi- 
derably. To calculate space-energy distributions of neutrons one needs to solve the systems 
of equations for subgroups by methods developed in thcrmalization problems 

But subgroups number is not found to be large. For evaluation of averaged parameters 
in Pj- and P 3 — approximations 2 and 3 subgroups correspondingly will do. 

CALCULATION OF AVERAGED CROSS-SECTIONS 
Cross-sections averaged in different ways were shown above to be necessary for determina- 
tion of mean group-parameters (15). These values can be directly derived from experimental 
ata on cross-sections dependence on energy in low energy range only were details of cross- 
seetton structure are available. The energy regions where cross-sections are low anddifficult 
to measure contribute mainly to values of these parameters. 

3 5 7 
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The problem arising in this connection is to reconstruct the curve of cross-section 
energy dependence on the base of available experimental data measured with poor resolution. 
Some theoretical concepts of cross-section structure are to be applied in this case. 

C KOSS -SECTIONS STRUCTURE IN THE REGION OF 
ISOLATED RESONANCES 

The exact solution of the problem mentioned is possible only in the region of isolated 
resonances. Energy dependence of cross-sections in the vicinity of isolated level is given 
by Breit-Wigner formulae': 

0 = [ o Q /(1 + x2)](€os2(/>1 - xSin 2</q) + = o r + a p (17a) 

a x = [ Oo /(i + x 2 )] iyr (17b) 

P, (g)IB" n + 2(13 - A -xl)/ 2(21 + l)(l + x2)l = Md ( ,V L l^)(/ L o„) 

L L ! L (17c) 

where a 0 = 4mA 2 [2J + D/2(21 + 1)1 r n /r, x = 2(E - E 0 )/r 

r n andr are neutron and total widths of resonance, T x = is reaction width, E 0 - resonance 
energy, X - neutron wave-length, 4>\ potential phase-shirt corresponding to spin and parity of 
level in question; a r - ’’reaction cross-section” determined in formal way. In (17b) B ^ 
describes the contribution of potential scattering to L th hnnnonic of scattering cross-section. 

B nn = 2 (2) 1 +l )(21 + l)(l 1 010L0) 2 Sin^ x Sin Cos( ^ - d-j) ( 1BaJ 

Coefficients B^, Al, II are determined by resonance parameters and phase-shifts of 
potential scattering interferring with resonance: % 

b l = z 2(UUiiL)(r2 /r 2 ) Cos2(d>ji — ^> " 

A . . 2(2J +1) S (211 + 1) (UOlO)LO)2(r n /T)Sin^ 1 Sin(2^^ 1 i) (18c) 
l L = 2(2J + 1)2 ( 2 t 1 + l)(l 1 010!1.0)2(r l )/r)Sin^ ll .Cos(2<Al-^ 1 i) U 8a > 

0 lOlO LO), zdjlll jL) are the Clebsh-Gordan and Blntt-Bidenhom coefficients tabulated 

in ref. 1181. , , 

The above formulae are valid in tliecase when the partial neutron l-wavc only contributes 

to resonance formation. The further assumption is that nuclei do not move. 

When thermal movement of nuclei is taken into account cross-sections formulae 17(a,b) 
will have a form [5,18]: 

a- a +a >]rrfW\ J co (^x /2 + i ^ /2) ixp(-2i^)] (19a) 

p o Y e i 

a = <7 V rr /4(r /DR [o> (£x/2 + if/2)] * <7 (T /D'H/.x) (19b) 

A U a L O A 

where £ = T/A, A = 2 >/k TE /A Doppler-width of level; - is the probability integral of 
complex argument, tabulated in details (see ref. 19). 

Equation (17) may be used instead of (19) only when £»l condition (2) holds. Doppler- 
broadening effects considerably on resonance's form in opposite case. 

Doppler-broadening of resonances is not of major interest for consideration of scattering 
anisotropy, because in the cases when this effect is appreciable scattering resonances with 
high orbital moments can not be in general treated as isolated. 
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AVERAGING OF CROSS-SE( iTJONS IN THE REGION OF 
RESOLVED ISOLATE!) RESONANCES 

In this energy range cross-sections (Kg 16;, ^m,g,:d over large number of resonances 
may be expressed in terms of mean resonance parameters and evaluated. 

The method for evaluation of average cross scions needed for calculations in diffusion- 
approximation has been discussed in detail in refs (12] and [151. We represent here more 
general expressions valid for calculations of cross-sections needed in higher approximations 
De rivation procedure is identical with one applied in references cited and is omitted here. 

<M<r + <,)"> = [!/ (a + a)l 5 C " I 

1 m=0 “ r - ra (20a) 

<0 ^L /( ~ a +<') > = (, p /l j p> l/(rr + o) n '> + J Ln (20b) 

‘r.m “ > CT r /(° +<7 ) m 1 9 1 n ) 


J Ln - >CT si4/ (o+a)n (21b) 

Define also = > ° x /(v + <*) > 

Substitution of Hq (17) or (19) into (21) results in 

t 1 v- - m + ,°° v Cos 2 ~xSin m 

r,m = IT" k 2E^k Vm ^ k ' “k’^kH {, )[ 1 d * (22a) 

k ok x 2 -xct k Cos2<6 k +^. Cos 2<£k + 1 


n. -x/s 2 lki + x 2)dx 


In particular 


Ln x ** n 17 . p a k J o n 

Au “‘-ok'Mik (x 2 - Xa k Cos <£ k +a k Cos<£ k +l) 


( 1 v gr k 71 (4- “k- ^k^C^^k-^k^Sin 2 ^^ ^ ^ 

r ' ! = Au k k 2K ok ak [d ~ a k Sin 2 <^)(l + a k Cos 2 <jJ k )]^ ^ 

(22c) 

Au k 4k ok [(1 _ „ k si n 2 q + ak (; 0 s 2 ^ k ) 3 2 

, u _ 2,72 2 £k^k yLk-/»LkT Si " 2 » k ( 

P k Eok y 1 (1 -a k Sin 2 <^ k )(l + a k Cos 2 (ji k ) 

o a \ 

r 2 „ C k A 2 3yLk( 2 + «kCos2<i) -^ L k d LkSin2<A k (4 + 3a k Cos2^ kT Sin 2 ^ 

; r-77— wr— = 2 (2RI 


1 ^ ff ' r xk 


2o 2 Au k E ok [(1 - a k Sin 2 0 k ) (I + n k Cos 2 

I 1 ^ Irr xk a c k' ^ 

A k 2L 0 k [(1 — a k Sin 2 <^ k ) (1 -a k Cos 2 <£ k )]'^ 
The sum is taken over all levels ”k” in interval of averaging Au 


uTT “k 


A c + b L „ v II 

“k = a ° k/a P + ° rL = f 2(21 + 1) ’ Pl= ] W7l) 
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The functions jjj and in absence of Doppler-broadening are equal to unity. 

These functions were calculated elsewhere [4, 5, 9] form = 1, 2 in assumption that there 
in no interference between potential and resonance scattering (<£ = 0). In this case ^ = ,,, 
There are no principal difficulties in calculations of these functions when <f> ^ 0. 

REGION OF UNRESOLVED ISOLATED RESONANCES 


Mean resonance parameters may be used for calculations of the integral characteristics (21) 
m the energy range where resonances may be treated as isolated but are not resolved 
experimentally. 

The cross-sections at fixed energy are the sums of partial cross-sectionsover possible 
combinations of spin and parity*: 


I =S1 

r > m „ 




■■ 

where v = J, „ . Values \ v , J v . , I* may be expressed as follows: 
r.m* L,m* x 

jv jy , — — — +°° Cos20 — xSin20 

r,m=^ Dl % ’'n. l( fo'“v^y»ofo'V/ l 


(25) 


~ x^-xajSin20 i< + a i/ Cos2s4 (/ +l 


l m dx 


(26) 


- 2 d) n-1 dx 


J- ( TL^-^Lo^ +1 

Ln n^i V J — I 

°° (x^ — xa^Sin2<^+a^Cos2<£ h- l) n 


Da" 

P 


\ v -Z — - 


IvV 


where 


X 2 ^ [(l-a Sin 2 tf )(1 + I Cos 2 <f> )]^ 

V T v V V 

a =4*A 2 [(2J + l)/2(2I + l))r /r („ +o+l o 1 

v ni/ n in I;-*- 


(27) 

(28) 
(29) 


* i/ p ^ 1=() p,u 

The sum in latter expression's denominator reflects overlapping of resonances with dif- 
ferent spins and parities. 

D 0 - is average level spacing for spiln J and parity n. All the parameters are evaluated at the 
mean energy of interval. Functions x^ f x m> correspond to the effects of fluctuations of 
partial width. These fluctuations are described by Porter- 11 mas distribution [20] and functions 
mentioned may be calculated numerically. The results of such calculations for some of these 
functions for may be found in refs [11, 21]. Another aspects of fluctuations effect on average 
values of cross-sections will be discussed below. 


REGION OF RESONANCES PARTIALY OVERLAPPED BECAUSE 
OF DOPPLER-BROADENING 

The calculations of averaged cross-section in the region of partialy overlapped interfering 
levels were impossible up to date because of the lack of satisfactory theoretical representa- 
tion of cross-section structure in this energy range. It is possible to estimate the selfindicati- 
on effects only when overlapping of resonances begins because of Doppler-broadening (A D , 
but r« D). 

The method for calculation of averaged cross-scctions for this case is developed. In this 
method averaged cross-sections are represented as difference between real average cross- 
section and some correction which takes into account sellshielding effects and has a dispertion 
form: 
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<o x /o'>/<\/o'> = cy x ~(l/a)(CT^a - u^a ) 

<1/ cr n *> => l/<T n + l/2n(n + l)(o^ -a?)o l] ^ (30) 

Supposing Doppler-width to be large and level spacing to be constant we calculate this 
correction 

W - <%o = ^a Xv a pi ,[2lxp(-2rr 2 /b^) + (b l/ / v /2rr)'^ /S x ] 

v x 

a 2 _ a 2 « £ 0 2[21xp(-2rr2/b2 ) + (b v /\]2v) (f> v ] (31) 

V V 

where v - J, rr, b u = D^/A, functions $.* ^ ^oc [10, 12, 21] reflect the effects of reduced 
width distribution on the cross-section value. 

THE EFFECTS OF RESONANCE-POTENTIAL INTERFERENCE 
ON THE RESONANCE SELFSHIELDING OF CROSS-SECTIONS 

The interference between resonance and potential scattering, results in sharp change of 

the form of resonance line. Cross-section' s maximum is not placed at resonance energy and value 

of this shift is E - = (r/2)tg^ h maximum cross-section o = o + o . . where 

inux o i max o min 

2 

"min = °p ~ a 0 Sin *| is rea ched at energy E min = E q - (r/2) ctg^,. 

For even-even nuclei which after absorption of s-neutron form compound-states with the 
single possible spin-value J = 1/2 the minimum of cross-section equals to.^p « (1 -_IL). This 

value may be very small, when the widths of inelastic processes are small. For example Fe^ 
cross-section in the interference minimum of 29 kev resonance equals approximately 0.2 barn 
and is almost intirely due to the impurities of other isotopes of Fe. 

The existence of this minima is evident by of major importance for determination of dif- 
fusion coefficient. But the value of effective resonance absorption integral is subject to the 
same effect in some cases. (Fig. 1). 

But interference between potential and resonance scattering may be especially important 
for calculation of the dependence of effective resonance integral on temperature. (Table 1). 

The same effect takes place for other resonances of with high energy and large 
neutron width. 

The effect in question may be considerable in the region of partially overlapped levels 
too (see Fig. 2). Because of interference the condition T>D does not necesserily imply that 
there is no resonance structure of cross-sections. Cross-sections become smooth only when 
large number of decay channel of compound nucleus are open. 

EFFECTS OF NONUNIFORM COLLISION-DENSITY 

As was noted above, Plackek oscillation of collision density does not effect strongly on 
averaged cross-section in the case of narrow resonances. However these fluctuations are to be 
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taken into account for resonances with widths which ure comparable with energy loss on 
the nuclei of the element considered but are essentially less than energy loss due to the 
scattering on nuclei of other mixture — components, To find energy dependence of collision 
density one must solve integral equation, 

The temperature dependence of absorption resonance integral (barns) for 663eV - level 
of U 238 0 = 0. 

Table 1 


L 

Temperature K° 

4 . . . .. 

0 

300 

900 

- 0 

0.011 

0.012 

0.014 

,/o 

0*027 

0.024 

0.021 


Effects of resonance selfshielding on the average neutron scattering characteristics of the 
ox yg en « n, _ L , _ 



No selfslue Iding 



V 

££s 

Medium 




V 

0.2-0.4 

3.80 

0.134 

-0.14 0 

0.4-0 .8 

5.60 

0.086 

0.23 0 

0.8-1 .4 

4.30 

0.110 

0.08 0 

1.4-2 .5 

1 

1.75 

0.100 

0.12 0 


Table 


Selfshielding coefficients 


<V2f 0t > 


(g— gi/2?ot^s 



V (u) - 1 f w (u 1 ) [ 1 Mu 1 - u)/2(u)] du (32) 

iu-rj i 

In some cases effects of nonuniform collision density may lead to considerable change of 
results of averaging. 

As an example we can take the effective resonance absorption integral as a function of 
dilution for 29keV level of Fe^. It is demonstrated in Fig.l. The change of this function relative 
to the case of reniform collision density is determined by arising of collision density peak in 
the region of interference minimum where ratio l c /2 increases sharply. Another example is 
energy dependence of collision density near the first resonance of at 6,7 eV, calculated 
in [11) by Monte- Carlo method. 
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EFFECTS OF RESONANCE SELFSHIELDING ON THE AVERAGED 
PARAMETERS OF SCATTERING ANISOTROPY 

To apply Eq. (22) to calculation of averaged values of scattering anisotropy parameters 
one needs to known potential phase-shifts and spins of resonance levels formed after absorp- 
tion of the neutrons with high orbital momenta. The poor information on this question is 
available up to date. Relatively complete and reliable data exist only for first levels of light 
nuclei. The calculations carriedout for these nuclei show that resonance selfshielding can 
effect considerably on average angular distributions and consequently on the neutron diffusion. 
Averaged characteristics of scattering cross-section on oxygen are represented as an example 
in the Table 2. 

The averaging procedure was carried out for group integrals presented in [16] on the bases 
of angular distributions measured by Langsdorf et al [22]. 

The results of calculations may be applied not only to pure oxigen but also to the oxigen 
nuclei included in UO 2 and II 2 O. The fluctuations of collision-density were neglected. 

The selfshielding factors of fission and radiative capture cross-sections of Pu^39 


Table III 


E 


p u 239 



ft 



Sc 


H 

x 0 

f r\ 

10 3 

10 2 

10 

0 

10? 

10 2 

10 

0 



300 

1.00 

0.99 

0.97 

0.94 

1.00 

1.00 

0.98 

EMI 

21.5—46.5 kev 

10 

900 

1.00 

1.00 

0.98 

0.97 

1.00 

1.00 

0.99 

Wm 



2100 

1.00 

1.00 

0.99 

0.98 

1.00 

1.00 

0.99 

0.99 




1.00 

0.98 

0.93 

0.88 

1.00 

0.99 

0.95 

EH 

10.0—21 .5 kev 



1.00 

0.99 

0.97 

093 

1.00 

1.00 

0.98 

IH 




1.00 

0.99 

0.98 

0.95 

1.00 

1.00 

0.99 

0.97 




1.00 

0.96 

mm 

0.80 

1.00 

0.96 

0.87 

ESI 

4.65+10.0 kev 

12 


1.00 

0.98 

wBm 

0.87 

1.00 

0.98 

0.94 

n 



2100 

1.00 

0.99 

EH 

0.91 

1.00 

0.99 

0.96 

0.92 



300 

ESI 

0.92 

ESI 

0.69 


0.92 

0.77 

0.67 

2.15-4.65 kev 

13 

900 

eh 

0.96 

EH 

0.79 


0.96 

0.87 

0.78 



2100 

LOO 

0.98 

EH 

0.88 

1.00 

0.98 

0.91 

0.87 



BEB 

ESI 

0.88 

EH 

0.56 

0.98 

IH 

0.62 


1.00—2.15 kev 


■ 

EH 

0.92 

IH 

0.69 

0.99 

EH 

0.76 

El 




1.00 

0.95 

EH 

0.82 

1.00 


0.84 

0.79 



300 

m 

0.80 

0.53 

0.42 

mm 


0.46 

0.36 

465—1000 ev 

15 

900 

E H 

0.86 

0.64 

0.53 

EH 


0.58 

0.47 



2100 

EH 

0.90 

0.77 

0.71 

0.99 


0.72 

0.64 



300 

0.93 

0.70 

IS! 

■H 

0.91 


0.33 


215-465 ev 

16 

900 

0.96 

0.78 

EH 

iH 

0.94 


0.42 




2100 

0.98 

0.83 

EH 

■H 

0.96 

0.80 

0.53 

0.45 



300 

0.91 

0.63 

0.34 

0.25 

0.89 

0.58 

0.28 

0.20 

100-215 ev 

17 

900 

0.94 

0.69 

0.41 

0.32 

0.92 

0.64 

0.34 

0.26 



2100 

0.97 

0.75 

0.49 

0.39 

0.95 

0.70 

0.41 

0.32 
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E 

i 

■ 

Pu 239 



ft 



fc 


I 



_L0 2 



10 3 

10 2 


0 




m 

0.50 


0.34 

mrm 

KOI 

PfS 

o 

46.5-100 ev 

18 

900 | 

n 


m 

0.42 

in 

liJ 

Ml 

M 


The proper theoretical analysis of resonance effects is necessary for construction of 
multygroup constant system for fissionable nuclei on the basis of available experimental data. 

We consider the effects of fission and neutron widths fluctuations on various average 
characterizing capture to fission ratio of probabilities. There exist different experimental 
procedures for determination of such ratio. One can measure a) cross-sections of capture a 

and fission <rf; b) effective number of neutron yield per neutron absorbed vo t Y 

1 = — 

CT f+ °.y 

(u - the average number of neutrons emitted per fission); c) capture to fission ratio (relative 

measurements) a y / af. Immediate averaging of these values over large energy interval results in 
following mean characteristics: 


1)a = <yaf, 2)a =i//,_ 1 , Where = > 3 ) < a ,> = ( JV_) 

' ty + nf <7f 

The last ration of this kind is a = f y /{* To obtain this one it's necessary to analyze 
cross-sections structure in resonance-region. All these ratios are d i f f e r e n t. Among 
them a is of greatest importance in reactor calculations. To demonstrate quantitative difference 
between "alpha’s" we consider a system of isolated noninterfering levels of compound- 
nucleus with the same spin and parity. Statistical fluctuations of fission and neutron widths 
will be taken into account in the way developed by Bethe [293 Oleksa [31], Lane and Lynn [21] 
orter and Thomas [20]. The procedure of averaging over and fluctuations we designate by 
brackets. Then (supposing not to fluctuate) we get 


a = r y > 


r n . r f r n , 

r y +iM' n > 7 < r f +r y +r n > - a = a v = r y < ■ 


>, a 


r f+r n ’ “ <r f > 

Using x 2 - distribution functions with the number of degrees of freedom which 
fission width and unity for neutron width one can see that 


, (33) 


equals u for 


> a^> a > a ( 34 ) 

a decreases approaching a when energy increases even with r, and T being constant - only 
because of neutron width increasing, before neutron p-were becomes important 

In fact s-neutrons absorbed by nucleus with spin I form compound - states with momenta • 

J = I + 1/2 and statistical weights proportional to 2 J + I. Thus cross-sections and widths 
must be averaged separatly for each level set with corresponding values of average fission 
and reduced neutron width and in general with respective number of fission channels. Eq (30) 

W 
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are to be modified correspondingly for this case, but (31) continues to hold. Resonance inter- 
ference and overlapping do not change considerably the area under resonance curve and 
consequently do not affect on a, i.e. do not breuk inequality a *> a which is to hold. Calcula- 
tions for Pu239 taking the existence of two resonance subsets into account were made and 
confirmed this conclusion (see Fig. 3). Values of parameters used in these calculations are 
indicated in Fig. 3 caption. 

For Lj233 a =0.15 and a = 0.35 [31]. Calculated values are a - 0.27 and 0.34 for v - 2 
and 3 correspondingly. Average fission widths for resonance subsets were taken from ref. [32]. 

But for [j235 (34) begins to fail approximately at 7eV and higher: for 30 resonances above 
this energy a * 1.4 and a = 1.2 [33], [34], [39]. This contradiction can be avoided by supposi- 
tion that some part of weak and broad resonances with small rg and large Tf have been missed 
in experiment. Formulae estimated the consequences of such missing quantitatively were 
derived in [35]. The explane mentioned and some other anomalies in resonance structure of 
L)235 i n the energy range under consideration — the existance of background in fission cross- 
section 30 barn / y/E , sharp difference between experimental distribution of reduced neutron 
widths and one - channel x 2 - distribution for small F° and relatively small observed value 
of mean fission width for resonances II - 40 TfQ^o) = 29 mev (for first ten resonances this 
value rfQ_ 10 ) = 78 mev, see Fig. 4). Real a calculated with level missing taken into account 
equals 0.85. This value satisfies (34) and corresponds to a (0.85) = 1.15 at u = 3. This value 
of a is more ill line with experiment. 

If this assumption is true, real value of average level spacing and mean reduced neutron 
width are approximately 1.5-2 times less than observed values. 

Integrals in the expression for <a> diverge at 1.2 (33). That means in fact that with 
such broad distributions of fission widths very small values <a> have considerable probability 

and Tf - value is determined by chance but exceeds a, a and a greatly. The integrals are 

. % J 
finite for y >2 and <<r> 3a and 2 a for a = 3 and 4 correspondingly. 

But interference and overlapping of resonances effect on <a> (namely decrease it) 
considerably. The approximation of isolated levels is justified in the case of Pu 239 and for 
this nucleus really <a> 2.5 » a 0.4 as expected for u = 1. For U 233 <a> 0.35 w>a = 0.15. 
The interference effects and overlapping are strong in the case of U 235 and < a > is close 
to a for this nucleus. 

The effects discussed above may result in following fenomenon. In relative measurement 
of o y / <rf with energy resolution better than level spacing a y / a ( r y /T f of corresponding 
resonance at every energy point and oscillates strongly reflecting cross-section structure. 

When averaged this value will give <a>. Selector resolution falls with encreasing energy and 
when it becomes worse than level spacing, cross-sections would be automatically averaged 
both in numerator and denominator and <a> would "turn into a” - the mean value of quantity 
to be measured decreases rapidly without slightest physical reason. It means that results of 
a y /o { measurement carried out with poor resolution - with energy spread being approximately 
equal to group interval _ (but with good statistics) - are of great interest because they give 
immediately a for the group under consideration. 
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In 26-group constant system published in [16] calculations of selfshielding factors for 
fission and capture Jf, / c for Pl» 239 were carried out with mean parameters from [36]. 

Fission width was supposed to be described by three-channel distribution as in [37]. Mean- 
parameters measured in [33] p239 an d one channel distribution are evidently more in line 
with experiment. They have been used for calculations of / c and /(• with results shown in 
Table 3. Below IkeV resonances were treated as narrow arid overlapping c;' '^vels with the 
same J was neglected. F n andTf fluctuations were taken into account also in [37]. The 
dependence of fission and capture cross-sections on temperature was determined by using 
functions calculated in [9]. 

Selfshielding coefficients in KeV-region where resonances begin to overlap because of 
Doppler-broadening were calculated by method developed in 10. The Table 3 containing /f 
and / c shows that in this region selfshielding for of is stronger than for a y with resulting 
increasing of a. 

Selfshielding factors were interpolated for 12-H groups. No new calculation of /f and f c 
were carried out for 18-th group. 

Doppler-coefficients for fast oxide reactor with volume 2000 liters were calculated with 
fc an d /f fr° m Table 3. Results are very close to ones obtained with coefficients published 
in [16]. 


EXPERIMENTAL DETERMINATION OF THE AVERAGED 
CROSS-SECTIONS 

The average cross-section calculations by using one-level formulae are valid only in rather 
narrow region of neutron energies where the effects of resonance interference and contribution 
of high angular momentum resonances (for which the average parameters usually are poorly 
known) can be neglected. As for higher energies the calculation of the quantities depending 
on the structure of the cross-sections becomes unreliable. Thus the problem of experimental 
determination of such quantities arises. 

Earliear [7, 12, 13] the experimental method of obtaining the parameters of the total cross- 
section structure had been suggested. By analysing the transmission curve one can obtain 
such qualities as for example <l/S n > (n = 1, 2, 3....). The measurements of the other 
necessary characteristics for the description of the neutron propagation in medium may be 
obtained in the same way. The assumption is made that the collision density fluctuations 
do not affect the results significantly. 


METHOD 

Collimated neutron beam with energy spectrum /(E) passed through sample of variable 
thickness. By using the AE detector with efficiency f (E) and holding ” the good geometry’ 1 
conditions one was allowed to measured the transmission function T(t): 

T(t) = N(t)/N(0) =/f(EME)r ME)t dli// r(EWE)dE (35) 

At, AE 

Here N - nuclear density of the sample. 
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Let us introduce P (a) - the total cross-section probability distribution. Then transmis- 
sion T(t). 

T(t) = J«[E W1P(a)l“ NB * da/jHE(a)lPMda (36) 

Mere E (a) is the function inverse to a(E); integration expances through the ail values of‘< 7 . 

If the energy interval AE used in transmission experiment is comparatively small one ccruld 
choose thedetector with nearly constant efficiency in this interval. The very interesting 
case is one of the detector sensitive to only one partial process (or one combination of 
partial processes) only, taking place in thin sample to be explored. Then the energy 
sensitivity cf the detector is proportional to - the cross-section of the chosen process and 
thus correlates with the variations of the total cross-section. 

Suppose we have the set of such detectors including the usual detector with flat ef- 
ficiency. The transmission curves measured by such detectors are the next 

T(t) = /P(o-)I' NQt d < r (37) 

T x (t) = /a x [EWlPWr Nat (38) 

T sL (t) = Jffg[E (<r)]/ L [E(<7)lP(a)l' N< ' t da (39) 

Now we make a representation of the cross section P(a) po babi Li ty distribution P(o-) as 
a sum of the 8 — functions': 

m 

PW- 2 a jfita-a.) (40) 

This is equivalent to the sunstitution of the real cross-section curve by the histogram in the 
energy interval AE. 

Then Eqs (37) — (39) become 

T(t) 


. £ ad-N-it 

(41) 

i=l 


m -No-t 
= S b:l 1 

(42) 

i=l ' 


m -No:t 

0 = 1 Cjl 1 

(43) 


where b L = * a ci /7 c c i - • 'x (*s V are the reaction (tissi ° n ’ Cap ' Ure CtC) 
cross-section and n th harmonic s of the elastic cross-section averaged over the neutron beam 
spectrum, a x and (<r s f L ) - the same quantities averaged over the spectrum of those part of 
beam neutrons for which average total cross-section has a value a s . It is easy to notice that 
the abovementioned quantities are just the same ones which are {resented in subgroup method. 
The set of transmission curves is necessary for this method. By the simultaneous least 
squares treatment of these curves one could obtain representation in the form (Eqs (4)). In 
the same way ids possible to evaluate for any curve the maximum number of the exponents 
in the curve expansion (justified by experimental accuracy) as well as uncertainty matnxe. 

Average cross-sections could be measured with the same detectors by convendonal 
methods. The values of parameters Eq.15 relate to subgroup characteristics in the way: 

<l/(a + a) tt >=ta i [l/(<T i + <r)] n (44) 

3 57 1=1 
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in , v 

<<jy/(o + o)> — X bj[l/(os + cr)l (45) 

i 9] 1 

in 

<a s^L + cr ) U> = ^ C|[ 1 /(aj-fcr)] 11 (46) 

U3 

Also these values can de obtained just from the transmission curves by using sucli relations 
as for instance. 

<a s f L /(a + a)«> =(N n /n/) ^ ? 0 T s (47) 

which are useful for evaluation of the maximum Power (<l/(cr + a)> n ) corresppnding to 
reasonable accuracy. 

THE MEASUREMENTS OF THE TOTAL CROSS-SECTION 
STRUCTURE CHARACTERISTICS 

The discussed method depends on availability of the neutron sources with sufficient 
intensity as well as rather efficient detectors. Therefore until it had been used mostly for 
studying of the total cross-section structure. Transmission experiments were carried out at 
the Van-Uraaf accelerator by using 1 (p, n) reaction as a neutron source. Neutrons emitted in 
backward and forward direction had the energies varying from lOkeV up to 3Mev. Measurements 
of the total cross section characteristics in eV region for heavy nuclei are under development. 
Neutrons of pulsed fast reactor (L3R) [23] are used by the time-of-flight technique. Banks 
of proportional counters (enriched to 96 % 3^) were used as a detector. 

In addition to earlier reported [14] data the new results are presented in Tables (2). 

Hie time-of-flieht results relates to energy intervals used in [16]. 

MEASUREMENT OF THE EFFECT OF RESONANCE 

SELFSHIELDING ON ANGULAR DISTRIBUTION OF 
NEUTRONS SCATTERED BY Fe 

This measurement is yet the only one relating to partial cross-section self- indication 
carried out by using aforementioned method. Fe was choosen because of the most prominent 
resonance structure of the total cross-section in Mev region. On the other hand this element 
has a great practical interest in connection with its wide utilization as a constructive and shield- 
ing material. Fast neutron beam of fast B.R-5 reactor after scattering by thin Fe sample was 
detected by Th 222 fission chamber with effective threshold nearly 1.5 Mev. The spectrum of 
detected neutrons was very close to the corresponding part of fission spectrum. Transmission 
function as well as angular distributions of neutrons scattered by Fe sample after transmission 
through different thicknesses were measured. By the compatible analysis of transmission curves 
received in this way the detectable neutrons were divided into two subgroups with the next 
characteristics. 

The results of Fable HI can not be explained by smooth cross-section variations only. 

Although this effect undoubtedly takes place the subgroup characteristic difference is mainly 
caused by resonance cross-section structure. One could notice that anysotropy of neutrons 
elastic scattering for energies 3—4 Mev caused mostly by resonance scattering. 

The results allow to demonstrate the importance of the exact account of the resonance self 
shielding when the problem of neutron propogation through the thick layers is to be solved, 
ror example the asymptotic relaxation length corresponding to E n = 3 .5Mev is equal to 7.95 cm 
when calculated in P^-approximation [26]. The same quantity determined by subgroup method is 
found to 8.40 cm. 

757 - 18 — 
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Figoi. The dependence of the reso- 
nance absorption integral for 29 keV 
level of Fe-^on the dilution of natu- 
ral Fe by the substance with large 
energy logs in scattering and cross- 
sections F 

Dotted, cifrves^were calculated 
neglecting fluctuations of collision 
density. The results corresponding to 
these fluctuations taken into account 
are shown by solid curves. Interference 
of resonance and potential scattering 
was: 1 )neglected, 2 Considered 

ill 9 


Fige 2. CroBB-sootion trend in 
vioinlty of two olose resonances 
of even-even nuoleuB with potential 
scat te ring supposed to be determined 
by S-wave only. 
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Fig.3» Calculated and experimental 
values ofo6 for Pu239are compared. 


Average resonance parameters are: r 

= 40 MeV, f f=160 MeV, = 0.6MeV wn 
fp3tf38] • Experimental distribution cor- 
responds to 0 = 1 for both sets of 
levels £35J • Experimental values of t/ 
are in agreement with such assumption 
within the limits of statistical errors, 


Fig. 4. Experimental gistogram (3) of 
C, statistical distribution of reduced 
neutron widths for lj235 resonances 
from 11-th to 40-th are compared with 
one-channol S £ - distribution (I) and 
with distribution calculated in 35 
which takes level missing into account 


- 20 - 


Approved For Release 2009/08/17 : CIA-RDP88-00904R0001 001 00043-3 






